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Abstract 



We prove an asymptotic formula for the number of Eulerian circuits 
for graphs with strong mixing properties and with vertices having even 
degrees. The exact value is determined up to the multiplicative error 
0(n~^/^"'"^), where n is the number of vertices. 



1 Introduction 

Let G be a simple connected graph all of whose vertices have even degrees. 
An Eulerian circuit in G is a closed walk (see, for example, [2]) which uses 
every edge of G exactly once. Two Eulerian circuits are called equivalent if 
one is a cyclic permutation of the other. It is clear that the size of such an 
equivalence class equals the number of edges of graph G. Let EG{G) denote 
the number of equivalence classes of Eulerian circuits in G. 

The problem of counting the number of Eulerian circuits in an undirected 
simple graph (i.e. graph without loops and multiple edges) is complete for 
the class t^P, see [3] . Thus this problem is difficult in terms of the complexity 
theory. Moreover, it should be noted that in contrast to many other hard 
problems of counting on graphs (see, for example, pLj, |9j), even approximate 
and probabilistic polynomial algorithms for counting the number of Eulerian 
circuits are not known in the literature. 

As concerns the class of complete graphs Kn, the exact expression of the 
number of Eulerian circuits for odd n is unknown (it is clear that EC{K„,) = 
for even n) and only the asymptotic formula was obtained (see [ID]): as 
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for any fixed e > 0. 

In [H the approach of [10] was generahzed. In particular, the asymptotic 
behaviour of the number of Eulerian circuits was determined for graphs with 
large algebraic connectivity. This class of graphs we mean as the class of 
graphs having strong mixing properties and it can be also defined as the class 
of graphs having large Cheeger constant (isoperimetric number) or sufficiently 
large spectral gap from 1 for the second largest eigenvalue of the transition 
probability matrix of the random walk on the graph and large degree of each 
vertex. The equivalence of these definitions was proved, for example, in [B]. 

In the present work we continue studies of [B], [ID]- We prove an 
asymptotic formula for the number of Eulerian circuits of graphs having 
strong mixing properties. This result is presented in detail in Section 2 of 
the present work and was given for the first time, but without proof, in [6]. 

Actually, the estimation of the number of Eulerian circuits was reduced 
in [1] to estimating of an n-dimensional integral which is close to Gaussian- 
type. We partly repeat this reduction in Sections 3, 8 of the present paper. 
In addition, we develop an approach for estimating of integrals of such a type 
in Sections 4, 6, 7. We prove the main result in Section 5. 

We note also that in [5] the asymptotic behaviour of the number of Eule- 
rian orientations was determined for graphs having strong mixing properties. 
Apparently, proceeding from the results of ^ and the estimates of the present 
work, it is possible to prove the asymptotic formula for the number of Eule- 
rian orientations given in [6] . In a subsequent paper we plan to develop this 
approach. 

2 Main result 

Let G be an undirected simple graph with vertex set VG = {vi, V2, ■ ■ ■ ,Vn} 
and edge set EG. We define n x n matrix Q by 




-1, {vj,Vk}eEG, 

dj, j = k, 

0, otherwise. 



(2.1) 
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where n = \VG\ and dj denotes the degree of Vj G VG. The matrix Q = Q{G) 
is called the Laplacian matrix of the graph G. The eigenvalues Ai < A2 < 
. . . < A,i of the matrix Q are always non-negative real numbers and the 
number of zero eigenvalues of Q coincides with the number of connected 
components of G, in particular, Ai = 0. The eigenvalue A2 = X2{G) is 
called the algebraic connectivity of the graph G. In addition, the following 
inequalities hold: 

Ti 

2mind,- - n + 2 < A2 < mind,-. (2.2) 

j n-1 j ' 

For more information on the spectral properties of graphs see, for example, [7] 
and [Uj. 

According to Kirchhoff's Matrix- Tree-Theorem, see jS], we have that 

t(G') = lA2A3---A„ = detMn, (2.3) 
n 

where t{G) denotes the number of spanning trees of the graph G and Mn 
results from deleting the first row and the first column of Q. 
We call the graph G as 7-mixing graph, 7 > 0, if 

the algebraic connectivity A2 = A2(G) > 7|l^G|. (2.4) 

In addition, we recall that (see, |j6j): 

• The class of 7-mixing graphs can be also defined as the class of graphs 
having sufficiently large Cheeger constant (isoperimetric number) or 
spectral gap from 1 for the second largest eigenvalue of the transition 
probability matrix of the random walk on the graph and sufficiently 
large degree of each vertex. 

• Almost all graphs (asymptotically, in some probabilistic sense) are 7- 
mixing. 

Let 

W = Q-^ = {Q + J)-\ (2.5) 

where Q is the Laplacian matrix and J denotes the matrix with every entry 
1. Let d = («!, . . . , an) e M"" be defined by 

a, = W,,. (2.6) 
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Let /3 = Qd and 

^1 = nfi - Y E ^^^^^ ' (2-7) 

k=l \ j=k+l J 

where dj is the degree of Vj G VG. Let 

R{0=tTiA{OWA{i)W), (2.9) 



where tr(-) is the trace fucntion, A(^) denotes the diagonal matrix whose 

;o cc 
M" 1 

Kronecker deha. Let = R{e^^^) 



diagonal elements are equal to corresponding components of the vector Q^. 
Let e^'^) = (e^'^ ei'^^) G M" be defined by ef = 6jk, where 6jk is the 



k=l 
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Finally, let 



where 



n ^i'^' (2.11) 



^ 4(rf, + l)2 2(ci, + l)(4 + l) 4(4 + 1)2' ^^'^^^ 

The main result of the present work is the following theorem. 

Theorem 2.1. Let G be an undirected simple graph with n vertices vi,V2, ■ ■ ■ ,Vn 
having even degrees. Let G be a 'y -mixing graph for some 7 > 0. Then 



EG{G) = (1 + (5(G)) C1G2C3G4 (y - 1) ! 2' 



77' /-f I n — 1 JT, — 1 

— —-Jl — — 



(2.13) 

where Gi,G2,G3,G^ are defined according to ([22D, (I23D, fl^TTUD . fim]) . re- 
spectively, dj is the degree of vertex Vj, t{G) is the number of spanning trees 
of G and for any e > 

|(5(G)| <Gn-l/'+^ (2.14) 
where constant C > depends only on 7 and e. 
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Proof of Theorem 12. II is given in Section 5. This proof is based on results 
presented in Sections 3, 4. 

Remark 2.1. For the case of the complete graph we have that: 

A2(ir„) = n, EK,, = tiK^) = n^-\ Q{K^) = nl, (2.15) 
where I is indentity matrix of size n, 

/3 = Qa = =^ = Ca = 1, (2.16) 

r, = i?(e-<*=)) = ^tr(A(e-<*^))A(e-<^))) = (2.17) 



n 



Thus we get that as n — > oo 



71 ( n — 1 ) 
2 



{v,,Vk}eEK„ ^ ^ (2.19) 



ri(ri — 1) 



Combining fl2.13p - f l2.19p . we obtain that the result of Theorem 12.11 for the 
case of the complete graph is reduced to (11. ip . 



3 Reduction to the integral 

A directed tree with root f is a connected directed graph T such that v G VT 
has out- degree zero, and each other vertex has out-degree one. Thus, T is a 
tree which has each edge oriented towards v. 

Let G be a connected undirected simple graph with n vertices fi, f2, . . . , fn 
having even degrees. Note that for every spanning tree T of the graph G and 
any vertex Vr G VG there is only one orientation of the edges of T such that 
we obtain a directed tree with root Vr. We denote by % the set of directed 
trees with root Vr obtained in such a way. 



5 



We recall that (see Section 4 and formulas (4.6), (4.7) of [1]): 

" rl 

EC{G) = n (- - l)!2l^^l-"+V-"5, (3.1) 

where for any r G N, r < n, 

S= n cosA.fe^ W (l + itanA,fe) rff, (3.2) 

where Ajk = C,j — C,k and 

Un{p) = {(^1,6, . . . ,a) e : 101 < P for all j = l...n}. (3.3) 

We approach the integral by first estimating it in the region which is the 
asymptotically significant one. In what follows, we fix some small constant 
e > 0. Define 

Vo = {^eU„{n/2): j^, - ^1 ( modyr) < ^-^2+^, ^ = ^^±^1^}, (3.4) 

n 

and let 5*0 denote the contribution to S* of ,^ G Vq. 

Under assumptions of Theorem 12. 1^ we have that as — )■ oo 

S = [1 + {exp{-cn^'))) Sq (3.5) 

for some c > depending only on 7. For the proof of (13. 5p . see Theorem 6.3 

of m. 

The sum over % in the integrand of (13. 2p can be expressed as a determi- 
nant, according to the following theorem of [12], which is a generalization of 
aforementioned Kirchhoff 's Matrix- Tree- Theorem: 

Theorem 3.1. Let Wjk < j, k < n, j k) be arbitrary. Define the n x n 
matrix A by 

A,, = ( , (3.6) 

the sum being over 1 < r < n with r ^ j. For any r with 1 < r < n, let 
Mr denote the principal minor of A formed by removing row r and column 
r. Then 

detMr = J2 n ^i'^' (3.7) 
where the sum is over all directed trees T with VT = {vi, V2, ■ ■ ■ , f„} and root 

Vr- 
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Using formulas fl3.2p and (13. 5p , Theorem 13.11 and the Taylor series ex- 
pansion of cos Ajfc and tan A^fcin the region Vq, one can obtain the following 
proposition: 

Proposition 3.1. Let the assumptions of Theorem \2. 1\ hold. Then as n ^ oo 
5o= (l + 0(n-i/2+6e^)2-i/V/2n-Metg Int, (3.8) 

Int= j expfzf^Qa-^fQf-^ ^% + \R{i)\dl {?,.<d) 

C/^(„-l/2+e) V {v,,vu}<^EG J 

where Q, a. and R{C) ai^^ the same as in (12.51) . (12. 6p and (12.90 . respectively. 

We prove in detail Proposition 13.11 in Section 8. Actually, this proof was 
implicitly given in |3] (see Lemma 5.3 of [1]). 

Thus, we get that to prove Theorem 12.11 it remains only to estimate the 
integral Int of (13. 8p . 



4 Asymptotic estimates of integrals 

We fix constants a,b,e > 0. In this section we use notation f = 0{g) meaning 
that I/I < cl^fl for some c > depending only on a, b and e. 
Let p > 1 be a real number and x G M". Let 

^ix,rj . (4.1) 



For p = oo we have the maximum norm 

||x|| = max |a;j|. (4.2) 

j 

The matrix norm corresponding to the p-norm for vectors is 

\\M, = ^^P\^- (4.3) 
x^o Flip 

One can show that for symmetric matrix A and p > I 

\\A\\>\\A\\,. (4.4) 



For invertible matrices define the condition number 

f,,{A) = \\Al-\\A-\>\\AA-\ = l. (4.5) 
Let / be identity n x n matrix and A = I + X he such a matrix that: 



A is positive definite symmetric matrix, 
\X,k\ < a/n, Xjj = 0, 11^-^2 < b. 



(4.6) 



Note that 

n 

p-ia' < Plh < WMoo = Pill = max^|A,-fc| = 0(1). (4.7) 



We recall that (see Lemma 3.2 of [5] ), under assumptions (14.61) . 

^^oo{A)=fIl{A) = 0{^i2{A)). (4.8) 

Using (14. 4p . (14. 6p . (14. 7p and (14.80 . we obtain the following lemma: 
Lemma 4.1. Let A satisfy ( [//.6] j. Then 

\\A-'\\^ = \\A-% = 0{1), (4.9) 

|X;,| = 0(n-^), (4.10) 

where 

X' = A'^ - I = A-\I - A) = -A-^X. (4.11) 
We use the following notation: 

< 9 >F,n= j g0)e^^Ue, (4.12) 

For r > let 

< 9 >F,r = < 9 >F,U„{rn^), (4-13) 

where 

f/„(p) = {(0i,e2,...,^n)GK" : 1^,1 <pforallj = l...n}. (4.14) 
We use functions of the following type: 

F{e) = -(Fao + H{e), (4.15) 
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where A satisfy f l4.6p . We consider the following assumptions on function H: 



H{9) < ci- 



n 



dH{9) 



dO 



< C2 



+ 



n 



For the case when i7 = we use notations: 

< g >n=< 9 >F,n, < g >r=< 9 >F,r, < 9 >=< 9 >- 



(4.16) 
(4.17) 

(4.18) 



Proposition 4.1. Let i7 C M" 6e such that Un{rirf) C Vt C Un{r2'rf) for 
some r2 > ri > 0. Let A satisfy ( [^.6p and assumptions ( [y^.i6p , ( [^^.i?] ) hold 
for some Ci, C2 > 0. Then 



< 1 >F,Q= O (< 1 >) , 



< Ol >F,n= - < 1 >F,n +0{n 

3 

< Ol >F,n= - < 1 >F,n +Oin 



-l+4e\ 



< 1 >, 



< 1 > 



and, for k ^ I: 



< OkOi >F,n= 0{n 

< OkOf >F,Q= 0{n 



<1>, 
< 1 >, 



< OkOf >F,n= < 1 >F,n +0{n 



-l+7e\ 



< 1 >, 



(4.19) 
(4.20) 
(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 



where F is defined by (JjJ^ (^''^d C3 = c^iri, r2,Ci,C2,a,b,e) > 0. 

In addition, for any vector p = {pi,P2, ■ ■ - Pn) ^ 1^"; IIpIIoo = 0{n^^/'^), 



I 
2 



J2 PAA~'),k <e''^P-'P'^'>^ > 



(4.26) 



j¥=k,j<n 



+ 0{n 



< 1 >, 



where {A ^)jk denotes {j,k)-th element of the matrix A ^. 
Proof of Proposition 14. II is given in Section 6. 
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5 Proof of Theorem 



2.1 



The Laplacian matrix Q of the graph G has the eigenvector [1, 1, ... , 1]-^, 
corresponding to the eigenvalue Ai = 0. Let Q = Q + J, where J denotes 
the matrix with every entry 1. Note that Q and Q have the same set of 
eigenvectors and eigenvalues, except for the eigenvalue corresponding to the 
eigenvector [1,1,..., 1]"^, which equals for Q and n for Q. Using (12. 3p . we 
find that 

t{G) = -\,\s---\n = ^. (5.1) 
Using (14. 4p . we get that 

n 

An = IIQII2 < IIQII2 < IIQIIi = max^|4-fc| = n. (5.2) 
Then, we find that 



\Q-% = ^<—. (5.3) 
A2 7n 



Using (12. 2p . we get that 



Ti — 1 

n-l>dj>X2—^>-f{n-l). (5.4) 



Consider the integral of Proposition 13.11 



Int = 

C/n(n-l/2+e) 



{vj,v^:}eEG 



(5.5) 

p = Qd, 

where A(^) denotes the diagonal matrix whose diagonal elements are equal 
to corresponding components of the vector and a is the vector composed 
of the diagonal elements of Q~^. 

Let define m = iUO),UO), • • • UO)) by 



9k = v/(4 + l)/2efc. (5.6) 
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Then we can rewrite f lS.Sp in notations of Section 4: 



where p= {pi,p2, ■ ■ -Pn), 



Pk = , (5.8) 



F{9) = -(FAe + H{e), (5.9) 



a'A9 = ^^(efqm, = -======Q,k, (5.io) 

^ A/(dj + 1)(4 + 1) 



= A', + «, (5.11) 

{vj,Vk}eEG 

n = {eeW: ^(9) e f/„(r2-^/2+^)}. (5.12) 
We aim to reduce, using Proposition 14.11 expression Int of (15. 7p to 



< 1 >= / 6-^"^^^^-= ^= = n V + l)/2- (5-13) 



Our argument is as follows: first we have to verify that all assumptions of 
Proposition 14.11 hold, then we will gradually get rid of the oscillating term 
e*^^, quadratic term ^'^^^^'> and the residual term — ^ ^ A^^. 

Further, we always use notation / = 0{g) meaning that |/| < c\g\ for 
some constant c > depending only on 7 and e. 

5.1 Assumptions of Proposition 14.11 

Combining ([53]), dEl, dESD, (15:T0|) and fl5J2|) . we get that 

A satisfy dH]) and f/„(rin^) C C f/„(r2n^) (5.14) 
for some a, b,ri,r2 > depending only on 7. 
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Let e^'^^ = {ef'\ . . . , Cn^) G M" be defined by ej*^'' = 6jk, where Sjk is the 
Kronecker delta. Due to the linearity of A(^) and tr(-), we find that 

Rim) = aofRm = ^se, (5.15) 

where 

i?,fe = tr(A(e-<^))Q-iA(e-<^))g-i), 

g ^ Rjk (5.16) 

v/(ci, + l)/2v/(4 + l)/2' 
We use the following inequalities for n x n matrices X, Y: 



|tr(XF)| < \\X\\hs\\Y\\hs, 
\\XY\\Hs<\\X\\Hs\\Y^h, 
where || ■ \\hs denotes the Hilbert-Schmidt norm, 



(5.17) 



1X1 



HS 



\ 



EEi-'f^*!'- (5-18) 

3=1 k=l 



Combining (Q, (jEl, fl57[5|) -( l5Tr|) . we find that 

eT^6 < \mi)Q-'\\Hs\m2)Q-"'\\Hs < 

< \\Am\Hs\m2)\\Hs\\Q-'\\i = mihimuQ-x < (5-19) 

<IIQII^IIQ'^ll^llfl||2||6l|2 = 0(l)||fl|h||6l|2. 

Using (1521), CTl . f l5lTD and f lCTl) . we obtain that 

i7(e"') < Rm) = aofRm = ^(i)!!^^!!^ = 

n 

Let {Q~^ A{^)Q~^) denote the (fc, /c)-th element of the matrix Q~^A{^)Q~^ . 
For any 1 < < n, we have that 

^ = 2tr (^Q^'AmA = 2tr (Ai^'WmQ-') = 
= 2dkiQ-'A{0Q-')kk + 2tr (Ag-iA(e)g-^^ 
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where A is the diagonal matrix with the diagonal elements I Ajj I = |A(e^'^))jj| < 
1 in the case of j ^ k and \Akk\ = 0. In particular, we have that 

IIAII2 < 1. (5.22) 

Since A(^) is diagonal matrix, we get that 

\d,iQ-'AiOQ-')kk\ = \dkAiOkk\\\iQ-%\\l < dk\m)\\2\\iQ-%\\l (5.23) 
where {Q~^)k is the A;-th column of the matrix Q^^. Note that 

UQ- )k\\2 

We also note that 



" <||Q-'||2||g(Q'^).||2 = ||Q-^||2. 



llA(oii2 = iigeiioo<2n|ieiioo. 

Combining ([El, (Q, dESl-dElSD, we get that 

|tr (aQ-'A{OQ-') I < n\\AQ-^A{0Q-^\\2 < 



(5.24) 
(5.25) 



<^l|Ail2||A(ai|2||Q- 



-1||2 



(5.26) 



dk{Q''AiOQ^')kk = 0{imoo. 
Using (El, O, (EH, fICT]) and (E21 for all 1 < A; < n, we obtain that 



ma. 



2\\S9\ 



dnm) 



de 



= o{i)mv 

o{n-')\\ei 



(5.27) 



For any 1 < /c < n, we have that 



d 



— Y A^, 



< 



4 5^|A,,|3 = 0(n)||el 



3 

00 



(5.2^ 



Combining (15. 6p and (I5.28p . we find that 



d_ 

de 



{vj,vi}£EG 



(5.29) 
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Using fHTTU]) . the fact that 



and fl5.4l) . we get that 



iQ'')jk = ^ (5.30) 



bloc = sup = 0(n-V2)||ga|U = 0(^-^2), (5.3I) 

l<fc<n ^y{dj + l)/2 

where vector p is the same as in (15.81) . 

Putting together fISAOD . dEUD, flOOD . f lOTj) . ( K29\t and flOT]) . we get 
that all assumptions of Proposition 14. II with data (15.81) - (15.121) hold for some 
constants a, 6, ri, r2, Ci, C2 > depending only on 7. 

5.2 Oscillating term 

Let define p^''^ = {p^i\p^2\ ■ ■ ■ iPn^) ^ ^" follows: 

pf = [ °' ^ ^ ^' (5.32) 

■' y Pj^ otherwise. 

Let 

k 

Note that 

F(<^\e) = F{e). (5.34) 
Using (I5.3ip and the fact that 

J-E^X' = 2A = 0(n-^)||5||oo, (5.35) 

we find that all assumptions of Proposition 14 . 1 1 hold for the case when we take 
function F^^\9) instead of the function F(9), vector p^^-* instead of vector p 
and for some constants a, b, ri, r2, Ci, C2 > depending only on 7. 
Note that 

een ^ \\e\u = o(n^). (5.36) 
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Taking into account fl5.31l) and using the Taylor series expansion, we get that 
giPfcSfe = e-5P'^' + ipi^Ok + 0{n-^/^+^'), eeQ. (5.37) 
Using fl53T|) and flCTj) . f lCTj) for F^''^^) and p<'=-i), we find that 



fc=l \ j=A;+l 
n-1 



fc=i V j=fc+i 



(5.38) 



iPk < Oke ^ >F(fc-i),n= 

n 

j=k+i 

Combining (K20\i . f lOTjl and flSlBSj) . we obtain that 

1 - f E PM-')^k) <e^^^^^'=' >^(.),^ + (5.39) 

j=k+l J 

+ 0(n-3/2+5-) < 1 > . 
Using fHlOj) . dEHD, f l530|) and f l53T|) . we note also that 

p,Pfc(A-i),fc = /3,/3fc(g-^),fc = 0(n-'), J ^ /c, 

Using (15341) . ( lOOjl and ( 09|) for = 1, 2 . . . n, we get that 

< e^'^^^ >^,f,= < 1 >^(„),^ +0(n-i/2+5^) < 1 >, (5.41) 
where Ci is the same constant as in (12. 7p : 



(5.42) 



Taking into account (I5.3ip and using the Taylor series expansion, we get 
that 

e--yk'l = i-^plel + oin-'^'^), een. (5.43) 
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Combining fl5:i5]) and flOT]) for F^''\ we get that 

< 1 >FW,n= - \pI^ < 1 >F(*-i),Q +0(^-2+4^) < 1 > . 



Using fICT]) . dEMD and f lCTj) for A; = 1, 2 . . . , n, we find that 
< 1 >F(r.),n= C2 < 1 >F,n +0(72-1+4^) < 1 >, 
where C2 is the same constant as in (12. 8p : 



-32 



fc=i ^ ^ fc=i ^ 



2(4 + 1; 



(5.44) 



(5.45) 



(5.46) 



5.3 Quadratic term 

Let define 6'' = {O'l, . . . , 9'^) G W as follows: 



e 



0, 



for j < k, 



6j, otherwise. 



Let 



F\e) = -(FAe + H\e), 



(5.47) 



(5.48) 



{vj,Vk}eEG 

In absolutely similar way as given in Subsection 5.1, we find that all assump- 
tions of Proposition 14 . 1 1 hold for the case when we take F^{9) instead of F{9) 
for some constants a, 6, ri, r2, Ci, C2 > depending only on 7. 
Note that 

F^{e) = F{e), (5.49) 

Ri^ie')) = ^{e'fR^ie') = e'^se', (5.50) 

where matrices R, S are the same that in (I5.15p . (15.161) . 

Combining (15. 48 p . the Mean Value Theorem and (I5.27p . we get that 



dR{^{9*)) 



Ok 



2 

oo' 



(5.51) 
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— * —tj —tj 

where 9* lies on the segment between 9 and 9 . Using fl5.36p . fl5.5ip and 
the Taylor series expansion, we get that 



= 1 + -Skk9l+ J2 S,k9j9k + 0{n-^+''), 9eQ. 



j=k+i 

Using ([52D, dEIS]), (EH, we find that 



Sjk < 9j9k >pk-i^ 

j=k+i 



(5.52) 



(4 + l)/2 

Since IK27\) imply that H^Hoo = 0{n~^), using (SSHl), we get that 



< I sup < 9j9k >F''-\Q I 52 I'S'jfel = ^ 

fc<j<n ^C^i (5.54) 



Combining (ICTj) . (ICTj) . fl532|) and (153^1) . we obtain that 

< 1 >F>=-^,n= (l + \Skk^ < 1 >p.,f, +0(n-'+^^) < 1 > . (5.55) 

Using (I53S]), (ESSD and (E35D for A; = 1, 2 . . . , ra, we find that 

< 1 >F,n= < 1 >F«,n +0{n-^+'') < 1 >, (5.56) 
where C3 is the same constant as in fl2.10p : 

C3^n(i.is.).n0-^5^^TT))- <^-^^' 

5.4 Residual term 

For a subset G of EG we define 

Fe(0) = -iP'A0 + He(0), 
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In absolutely similar way as given in Subsection 5.1, we find that all assump- 
tions of Proposition 14. II hold for the case when we take Fq{6) instead of F{6) 
for some constants a, b,ri,r2,Ci,C2 > depending only on 7. 
Note that 

FEcie) = F'^ie), (5.59) 

~ (dj + ly " ^d,- + 1)3/2(4 + 1)1/2 + 1)(4 + 1)~ .ggQ. 

(4 + l)='/2(^^. + 1)1/2+ (4 + 1)2' 

Combining KT^ . (|4J4]), (g^SD, (El and (ICTj) . we get that 

4 _ /3 4 6 4 

< ,fc>Fe,n- I 4(rf. + i)2 + 4(^. + i)(4 + i) + 

^ (5.61) 

+ i(4TTFj<^>-e,.+0(.-3-)<i>- 

Using dOO]) and flS^GT]) . we get that 

< >p^^^= p^., < 1 >^^^^ +0(n-3+7-) < 1 >, (5.62) 

where 

" ^ ~ 4((i, + 1)2 ~ 2(d, + l)(4 + l) ~ 4(4 + 1)2' ^^'^^^ 

Note that 

1 - 1 < P,, < 1. (5.64) 

Using fl4.20p , f l5.62p , we can gradually remove all the edges from the residual 
term Heg and obtain that 

< 1 >F^a,^= n < 1 >^ +0{n'^+'') < 1 > . (5.65) 

Combining fl4.19p and (15. 65 p . we get that 

< 1 >Feg,^= < 1 > +0(n-i+^^) < 1 >, (5.66) 
where C4 is the same constant as in (12. lip . 
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Combining fICT]) . ( EM\f . (KM\f . flCT]) and (EM\f . we find that 

< e'^^ >F,n= C1C2C3C4 < 1 > +0(n-i/2+7s) < 1 > . (5.67) 

Putting together ([3ll]), ([33]), ([3SD, dEI]), dEZD. fl5J[3|) and fl537l) . we 
obtain (I2.13P and (12.141) (with the exponent 7e instead of e). 



6 Proof of Proposition 14.1 



In this section we use notation / = 0{g) meaning that |/| < c\g\ for some 
c > depending only on ri, r2, Ci, C2, a, b and e. 
Let 

m = (Md), 02 . . . , Md)) = AO- (6.1) 
According to (14.61) . A = I+X, Xjj = 0, and so for some gi{0) = (71 (6^2, . . . , On) 

e^Ae = <p\0) + g^(e). (6.2) 

Using (14. 6p . (16.21) and estimating insignificant parts of Gaussian integral of 
the following type: 

'"(max{|x|,fci})'e-("-'=^)'rfx, (6.3) 
we find that for r > 0, s > 

<ll^li:o> = / \\er^e-^^'de = 

+00 +00 / +00 



y \\er^e-f^^'ydeAde2...den (6.4) 

-00 —00 \-oo / 

= (1 + O (exp(-C4n2^))) j W^e-^^'dl 

\4>i{0)\<rn'' 

where C4 = C4(r, e, s) > 0. Combining similar to (16. 4p expressions for 0i, 02, • • . 4>n, 
we get that 

/ \\9\\l^e-'^^'d9= (1 + (exp(-C5n2^))) < >, (6.5) 

1 1^(0) 1 1 00 <r'n= 

19 



where C5 = C5(r, e,s) > 0. Combining fl4.9p . fl6.ip and fl6.5p with s = 0, we 
obtain KT^ . 

Using f l4.16p . we find that 

I < 1 >F,n I < J \e^^'^\de< J e-^^'+'^^^'de = 0{< 1 >). (6.6) 

In order to prove fl4.2ip - fl4.25p we use the following two lemmas. The 
proofs of them are given in Section 7. 

Lemma 6.1. Let Vt dW^ he such that Un{Tin^) d C. Un{r2n'^) for some 
r2 > ri > 0. Let A satisfy lji4-(>\ ) ('•''^d assumptions ( [y^.igp and ( [y^.i7[ j hold for 
some constants Ci,C2 > 0. Let P = P{x) = 0{\x\'^) for some fixed s > 0. 
Then for any T{e) such that \T{e)\ < P{ ^ 



< Tie) >Mn\n= O (expi-cen^n) < 1 >, (6-7) 
and for any T{6) = T{9i, . . . , Ok_i, 9k+i, ■ ■ ■ , On) such that T{6) < P{ 



+0 (exp(-C6n^^)) < 1 >, 

< 4m9) >F,n= I < m >F,n +0(n-i+^^) < \fie)\ >p,o + 

+0 (exp(-C6n2^)) < 1 >, 

< M0)T{9) >F,n= 0(n-i+^^) < \f{9)\ >F,n + 

+0 (exp(-C6n2^)) < 1 >, 

< (PU9)f{9) >F,n= 0{n-'+^') < \f{9)\ >F^n + 

+0 (exp(-C6n2^)) < 1 >, 



(6.10) 



(6.11) 



where function F is defined by fi4.15\ ), vector (f){9) is defined by lid. 1\) and 
constant cq = c^^ri, r2, Ci, C2, a, b, e, P) > 0. 

Lemma 6.2. Let assumptions of Lemma \ 6. 1\ hold andsi, S2, ■ ■ ■ , s„ G NU{0}, 

s = si + . . . + s„ > 0. 
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Let 

Sfc = and |{j : ^ 0}| < 3. (6.13) 

Then 

< MO)M{m) >F,n= 0(sn-i+(^+^)^) < 1 > . (6.14) 
Using ( KTU\i . we find that 

Ok = cj)k + qI$, \\q\\oo = Oin-'). (6.15) 
Combining (16. 6p . (I6.15p . Lemma [6. II and Lemma \6.2\ we obtain that: 

<6k{9f >F,n=< {ql$f >F,n= 

= Oin^') ( ^ < 0,(e? >p,n + Yl I < ^nWjM >F,n I ) = (6.16) 

= {0{n-^) + 0{n-^+'")) < 1 >= 0{n-'+^') < 1 >, 

< Skiey >F^n =< {fJY >F,n= 0{n-') ^ < 0,(0^ >F,n + 

j 

+ 0{n-^) J2 I < ^P.M^jM' >F,n 1 + 
ii7^i2 

+0(n-^) 5^ |<0,,(e"')0,,(^"')0,3(^"')2>^,n|+ (6-17) 
+0(n-4) ^ I < <p^,{9)<P,M^,M<PjM >p,n I = 

= (0(n-3) + 0(n-2+4e) ^ 0(^-3+7.) ^ ^^^_2+7e)^ 

+ 0(n-^+^^)) < 1 >= 0(72-^+^^) < 1 >, 

where 

6kiO) = ek-Md)- (6.18) 
According to (14. 6p . we have that 

4(^"') = 4(^1, . . . , 4-1, 4+1, ...,0n) (6.19) 
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Using dSSD, dSlS]), dSlS]), (EIIZD> fICT]) and Lemma EH we obtain that: 

< <Pk{0)h{0) >F,n = 

= 0(n-i+^^) < \6km >p,n +0 {expi~cen'')) < 1 >= (6.20) 
= 0{n-^+^') < 1 >, 

< 0l >F,n=< {MO)+Sk{e)f >F,n= 

=< >F,n +0(n-i+5^) < 1 >= (6,21) 

= ^ < 1 >F,n +0(n-i+5^) < 1 > . 

< MO)Sk{0)^ >F,n = 

= 0{n-'+^') < \6km' >F,n +0 {expi-cen'')) < 1 >= (6.22) 
= 0{n-^+^') < 1 >, 

1 



(6.23) 



+0 (exp(-C6n^^)) < 1 >= 
= 0{n-^+^') < 1 >, 

< MOf^kie) >F,n = 

= 0(n-i+6^) < Me)\ >F,n +0 {exp{-c,n'')) < 1 >= (6.24) 
= 0(n-^+^^) < 1 >, 

< et >F,n=< {<Pk{0)+h{e)Y >F,n= 

=< MOY >F,n +0{n-'+'') < 1 >= (6.25) 
< 1 >^,n +0(r2-i+^^) < 1 > . 



\4 ^ I nf^-i+7t 

_ 3 
~ 4 

In a similar way as in fl6.16p . f l6.17p . using fl6.20p - f l6.24p . we find that 



< 6u{e)ei >F,n = O(n-i) < 0,(^>z >F,n= 

^ (6.26) 
= (0(n"^+2^) + 0(n-^+^^)) < 1 >= 0(n-^+^^) < 1 >, 
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j 

+ Oin-') J2 I < ^PnWjMdf >F,n I = (6.27) 
= (0{n'^+^') + 0{n-^+'')^ < 1 >= 0(71-^+'') < 1 >, 

< Sk{e)ef >F,n= 

= 0{n-') J2 < U¥f >F,n +0{n-') < ^6)6^ >^,^= (g 28) 

= (0(n-^+^^) + 0(71-'+"')) < 1 >= 0(71-'+'') < 1 > . 
Using (ESD, (1619|) . dn^SD-dn^Hl) and Lemma EH we obtain that: 

< MO)h{0)Of >F,n= 0(n-i+4^) < M9)e^\ >F,n + 

+0 (exp(-C6ri2^)) < 1 >= (6.29) 
= 0{n-^+^') < 1 >, 

< 0k9i >F,n=< (MO) + kmOi >p,n= 0{n~'+'^) < 1 >, (6.30) 

< OkOf >F,n=< (Md) + hmef >F,n= 0(n-i+^^) < 1 >, (6.31) 

< ele^ >F,n =< (MS) + h{o)yef >F,n 

=< MOfO^i >F,u +0(n-i+^^) < 1 >= 

= 1 < Md) + m? >F,n +0(n-i+^^) < 1 >= ^^-22) 

= ^ < 1 >F,n +0(n-i+^^) < 1 > . 
Since ||pl|oo = O^n"^^"^), using the fact that 

ehpM = 1 + 0(n-i+2.), een, (6.33) 

and (16. 6p . we get that 

< Oke'^P-'P"'^'' >F,n=< Oke'^P-'P^'^'^ >F',n +0(n-^+=^") < 1 >, (6.34) 

where F' = F — \p\0'j.. It is clear that F' satisfy all assumptions of Lemma 
16.11 For any = (p'i,P2, . . . ,Pn) G such that llp'Hoo = 0{n~^^'^), we have 
that 

^ip'A = 1 + + 0(n^i+2"), e~'P''^' = 1 + 0{n-^/^+'), OeVt. (6.35) 
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Using dSSD, O and (KWf with F' instead of F, f05]) . we get that 

< 0;(^)e^^^P"-'fJ^' >p,,f,= 0(n-i+^^) < 1 >, 
< MWdy'^"^''"'''' >F',n= 0(n-i+5^) < 1 >, 
< Mme''^''"'''^'' >F',n=< (ff{e)e'^^-'^'fi' >F',n + 

= ^ < e"^^-'^'^'^ >p>^n +0{n-^+'") < 1 >= 

= i < e^'^^^ >,>n +0{n-'/'^'-) < 1 > . 

Combining fl6.35p and fl6.36p . we find that 

< MO)e''^'' >F',n=< MO)e''^''-'''^'' >F',n + 

+ < <Pi{e)eie''^^-'P''^ n +0(n-i+3-) < 1 >= (6.37) 

= ^P'l < e'^^^ >F',n +0{n~'+'n < 1 > • 

Using (US]), ([61]) and flOT]) . we get that 

n 

1=1 

Combining fl6.34p and fl6.38p for ff = p — pk(^^\ we obtain fl4.26p . 



(6.38) 



7 Proofs of Lemma 16.11 and Lemma 16.2 



In this section we continue use notation / = 0{g) meaning that |/| < c\g\ 
for some c > depending only on ri, r2, Ci, C2, a, h and e. 
Let 

^"('^^ = (^^i,...,e,_i,0,4+i,...,e„). (7.1) 
Proof of Lemma \6.1\ Using ( 14. 9 p , ( 16. ip and (16. 5p , we find that 

I < T >^.\n I < J PmUe-'^^'de = O {expi-cen'^)) < 1 > . (7.2) 
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For simplicity, let = 1. Using fl7.2p . we get that 



p= 1,2,3,4. 

where we let H{e) = for ^ G M" Combining (HTTll and the Mean Value 
Theorem, we find that for 9 G Un{r2n'^) 

H{e) - H{e^^^) = 0{n~^+^'). (7.4) 

Using (16. 2p . we get that 

1 0?(^"*)r(02, . . . , ^„)e-^^^^+^(^^)rf^"'+ O (exp(-C6n^^)) < 1 >= 

f^-gi{e2,...,en)+H{0W) 



(7.5) 



—r2n^ —r2n 



-r2ri' 



p = 0, 1,2,3,4. 
Combining ( K9\\ and ( TTil) . we find that for 6^^'^ G Un{r2rf) 



r2n 



—T2n'^ 

+ 00 +00 



y ^^e-'^?^^^)^;^! + O [eM-cir/')) j e'^"^'^de,+ (7.6) 

— oo ~oo 

+ y" 0?e-^?(^^) (e^(^>^(^^''') - l) rf^i, p = 0,1, 2, 3, 4, 



|</.i{0)|<r3n'= 
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where C7 = Cj{r2, Ci, C2, a, 6, e) > 0, rs = r^{r2, Ci, C2, a, 6, e) > 0. 
For p = 2, 4, we have that: 



+00 +00 



2 



-00 —00 
+00 +00 

(,4 



2/m .„ 3 / 

4 



-00 



\(pi{e)\<r3n^ 

+00 



-00 

+00 



0(n-i+"^) / e 



1) 



for G f/„(r2n^), p = 0,2,4. 



Combining Q-dH]), we obtain (EID and flOj) . 

For p = 1, 3, we have that: 



+00 

^2 



(7.7) 



(7J 



<Ple-^"^'^dei = (7.9) 
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|</)i(6l)|<r3nE 

0<<f>i{0)<r3n^ 



r-3n'=<flii(e)<0 



</.i>0 

+00 



for G t/„(r2n"), p = l,3. 
Combining (OD-dZS]), dZS]), (HH), fmU]) . we obtain ^Ml and (E3]) 

Proof of Lemma \6. 2[ Let T(^) satisfy 

\m\ = omu, 

^ip. = 0{sn-'-^)snp\n9)\, 9 e Q. 
Combining f l6.10p and (17. lip , we get that 

< (t>k0)T{e) >F,n= 

=< MOme^'^) >F,n + < M0)iTi9) - T{9^''^) >F,n= 

=< MOme^'^) >F,n +0{sn~'-') < sup {(PkOkTm >F,n= 

een 

= 0(sn-^+^^) < sup\T{9) \ >F,n +0 {exp{-cen^')) < 1 > 

een 

Using f l4.6p . we find that for 6* G f2 
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Combining f l6.6p andf l7.12p . we obtain f l6.14p 



8 Proof of Proposition 13.1 



In this section we use notation / = 0{g) meaning that |/| < c\g\ for some 
constant c > depending only on 7 and e. 

The following lemma will be applied to estimate the determinant of a 
matrix close to the identity matrix I. 

Lemma 8.1. Let \\\\ denote any matrix norm. Let X be such annxn matrix 
that \\X\\ < 1. Then for fixed m> 2 



/m-l /iNr+l \ 

det(J + X) = exp ^ ^ triX"-) + E„,{X) , (8.1) 



, r=l 



where tr{-) is the trace function and 

II '5^11™' 

|i^n.(X)|<-^^. (8.2) 

m 1 — \\X II 

The proof of Lemma 18.11 is based on estimation the trace of the matrix 
ln(/ + X), using the representation as a convergent series. Lemma (8.11 was 
also formulated and proved in [10]. 

For any 1 < r < n, we have that 

^0=/ n cosAjkJ2 n (l + ^tanA,fc) (8-3) 

where A^-fc = — and 

Vo = {^eUrX^/2): |e^.-^-|(mod7r)<n-V2+-,e=^i±lll±^}. (8.4) 

n 

Since the integrand is invariant under uniform translation of all the ^j's mod 
TT, we can fix ^ = and multiply it by the ratio of its range tt to the length 
n~"^/^ of the vector ^[1, 1, ... , 1]-^. Thus, for any 1 < r < n, we get that 

SQ = Tm^''^ I JJ cosAjfc^ JJ (1 + itanAjfc) dL, (8.5) 
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where L denotes the orthogonal complement to the vector [1,1,..., 1]"^. 
Let define n x n matrix B by 

— tanAjfc, for {vj,Vk} E EG, 



B,k={ ^ io^'^ = -^' (8.6) 



tanAj;, ioik = j 

:,Vi)(^EG 

otherwise . 

Using Theorem 13.11 with the matrix Q + iB, we get that 



n 



EE n (l + ^tanA,,) = ^M„ (8.7) 

r=l T&Tr (vj,Vk)£ET r=l 

where denotes the principal minor of A formed by removing row r and 
column r. Since the vector [1, 1, ... , 1]^ is the common eigenvector of the 
matrices Q and B, corresponding to the eigenvalue 0, we find that 

VM,^^^^^. (8,8) 

r=l 

where Q = Q + J and J denotes the matrix with every entry 1. Note that 
|A.fc| = l(0-e)-(a-e)l<2r^-l/'+^ ^ e Vo, (8.9) 

n 

1 15| |i = max = 0(^1/2+"), ieVo. (8.10) 

k=i 

Let $ = BQ-\ Using (jMl), (El, (lOnll and (18A0|) . we get that 

||*||i<||i?||i||Q"'||i = 0(72-1/2+^), ^eVo. (8.11) 

Using Lemma [8. II with the matrix i$, we find that 

det(J + i<l>) =exp + ^ + 0(72^^/2+^") j , ^ e Vq. (8.12) 

Let 

B = Bskew + Bdiag, (8.13) 

where Bskew is the skew-symmetric matrix and Bdiag is the diagonal matrix. 
Since Q is the symmetric matrix 

il{Bsken^Q-^) = 0. (8.14) 
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Using fl8.9p . we find that 

\B,,„. - AII2 = 0(n-i/2+3.)^ 



(8.15) 



where A denotes the diagonal matrix whose diagonal elements are equal to 
the components of the vector Q6. Combining f l5.3p and f l8.15p . we get that 



tr ((5,,,, - A)Q-i) I < n\\Bd,ag - M\2\\Q-% = 0(^-^/2+3^), 

Using (I8.14P and fl8.16p . we obtain that 

tr<l> = tT{Bd,aaQ-') = tr(Ag-i) + Oin-'/^+'') = 



(8.16) 



(8.17) 



where a is the vector composed of the diagonal elements of the matrix Q ^. 
Using the property of the trace function 

tr(XF) = tr(FX), (8.18) 

we have that 

tr<l>2 = tT{Bske.^Q-'f + tIiBd^agQ-'f + 2 tr (^Bsken^Q-' B^iagQ'') • (8.19) 

Since Bskew is the skew-symmetric matrix and Q^^B^iagQ^^ is the symmetric 
matrix, we find that 



tr (yBskewQ ^BdiagQ ^ j — 0. 

According to (15.170 . we have that 



(8.20) 



trX^ < 
|-^^||_f/5 ^ I l-^l IhsI |^"^| b- 



(8.21) 



Therefore we get that 



tr (-Bsfceii)Q 



-1\2 



^ \ \BskewQ Whs- 



(8.22) 
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Combining fl5.3p and f l8.9p . we obtain that 

mken^Q-'WHS < \\Q-%mkeJ\HS = 0(^-1/^+^), f ^ Vq. (8.23) 



Using ([53D, <B> and f lgJf^ . we get that 



tr [{Bdiag — A)(5 ^BdiagQ ^ 



< 



< n\\Q-'\\l\\{Bdiag - A)\\2\\Bdiag\\2 = 0(^-1+^^), ^ G ^/q, 

and 



(8.24) 



< 



^.25) 



tr l^iBaiag - A)Q-\Bdiag ' A)Q"'' 

< n||Q-l^||(5,,„, - A)||^ = 0(n-2+6^), f e K.- 

Combining (I8A8|) . (EMD and (K2^ . we obtain that 

tr{B,,,gQ-'f = ti{AQ-'f + 0(n-i+^^), f G K)- (8.26) 

Combining f lHTTOj) . ^Ml, (ESS), (E23]) and (E26]), we obtain that 

tr<l>2 = ti{AQ-'Y + 0(n-i+^^), f e K)- (8.27) 

Using flSTTj) and (18:271) in (18121) . we get that 

det(/ + «$) = exp (^t^Qa + ^^^^^ + 0(n-i/2+4^) j , ^ e K). (8.28) 

— * 

By Taylor's theorem we have that for ^ G Vq 
Y\ cos Ajk = 

(8.29) 



exp 



2 ^ •'''12 

(i,j,t;fc)G£;G (»;,,t)fe)G£;G 



Note also that 

E ^'. = e"*'<5f• (8.30) 
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Putting together ^ for r = 1, 2 . . . , n, IKm . ( K^f . IKMf . 

we obtain that 

5o = (1 + O (n-i/2+^^)) im-^/^detQ Int', (8.31) 

Int'= j expLc'Qd-^C'Q^-^ Yl ^'k + l^iOjdL, (8.32) 

where R{0 = tT{A{OQ-'A{OQ~')- 

Let -Pr(^) be the orthogonal projection ^ onto the space L, where L 
denotes the orthogonal complement to the vector [1, 1, . . . , 1]-^. Note that 



ii^^(Oiioo = lie - 1, • • • , 1]' Hoc < 2|ieiu, (8.33) 

where 

e=^i±^^±^. (8.34) 
n 

Thus 

Uniln-'/'+n C {( : Pr(a G f/„(n-V2+.)| (3.35) 
We also note that 

Q^=QPr{0- (8.36) 

Therefore the integrand of fl8.32p does not change under the substitution of 
vector e by vector Pr(e) and 

{vj,Vk}€EG / 

/+00 
J e-'^^^"dx = (8.37) 



exp I t 

LnVo 



P(e)eLnVo 



27r 

P(0GC/n{n-l/2 + £:) 



where 



m = -\i'Qi-To H + (8.38) 



12 

{vj,Vk}eEG 
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Using notations (EH]) - flHTTT]) . and formulas fOU]) . ([H35]), dSSZD, we get that 



where cy > depends only on 7 and e. 

Combining flHTSTjl and flHTSOD . we obtain flSTSj) . 
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